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Abstract 

We study the competing effects of gluon self-coupling and their interactions with 
quarks in a baryon, using the very simple setting of a hermitian 1 -matrix model with 
action tr A 4 — log det(z^ + A 2 ) . The logarithmic term comes from integrating out N 
quarks. The model is a caricature of 2d QCD coupled to adjoint scalars, which are 
the transversely polarized gluons in a dimensional reduction, v is a dimensionless 
ratio of quark mass to coupling constant. The model interpolates between gluons in 
the vacuum [y = oo) , gluons weakly coupled to heavy quarks (large v) and strongly 
coupled to light quarks in a baryon (u —* 0). It's solution in the large- N limit 
exhibits a phase transition from a weakly coupled 1-cut phase to a strongly coupled 
2-cut phase as v is decreased below z/ c = 0.27. Free energy and correlation functions 
are discontinuous in their third and second derivatives at u c . The transition to a 
two-cut phase forces eigenvalues of A away from zero, making glue-ring correlations 
grow as v is decreased. In particular, they are enhanced in a baryon compared to 
the vacuum. This investigation is motivated by a desire to understand why half the 
proton's momentum is contributed by gluons. 
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1 Introduction 



Photons carry a negligible amount of momentum in a moving atom. By contrast, it is experi- 
mentally found that gluons carry about half the momentum of a nucleon when observed at 
momentum transfers of order Q 2 ~ 1 (GeV) 2 . The growth of the gluon momentum contribution 
as Q 2 is increased is correctly predicted by perturbative QCD. However, the momentum frac- 
tion (xbj) dependence of gluons in a nucleon at any fixed value of Q 2 appears to be essentially 
non-perturbative. More generally, determining the 'emergent' bound-state structure of gluons in 
baryons from QCD remains an interesting and challenging open problem of theoretical physics. 
This problem in 3 or 4 dimensional QCD is still very hard to address analytically. However, 
there is a context where it can at least be finitely formulated. This is in the class of theories 
where 2d QCD is coupled to adjoint scalar fields. The principal examples are the reduction of 
QCD from 3 or 4 to 2 dimensions, where we assume that fields are independent of transverse 
coordinates . In this case, the longitudinal component of the gauge field contributes a linear 
potential between dynamical quarks q a {x) and the transverse components transform as adjoint 
scalar fields A%(x) , for instance in the light-cone gauge. Thus, for ./V colors a, b, we have a two- 
dimensional field theory of N x N hermitian matrix- valued and iV -component vector- valued 
fields with only a 'global' U (N) invariance. The gauge-invariant observables include glue-ring 
(closed-string) variables 



tr[A(x 1 )A(x 2 ) ■ ■ ■ A(x n )} 



(1) 



and meson (open-string) variables 



jI(x ) ah( X1 )a:i(x 2 )- 



(2) 



where the points Xi are null separated 2 . In 't Hooft's large- N limit, the expectation values of 
products of these variables factorize. So we may restrict attention to the single trace observables. 
Quarks are suppressed by one power of N in calculating vacuum expectation values of glue-ring 
observables in the large- N limit However, when the vacuum is replaced by a baryon state, 
interactions with quarks are just as important as gauge field self-interactions [B]. For definiteness, 
let us consider 2d QCD coupled to a single adjoint scalar field via the action (non-dynamical 
fields have been eliminated) 



S 



dt dx tr 



q'idtq + \m 2 q ] ^—q + d x Ad t A 
2 io T 



y / dtdxdy J fe a (x) 



y\ T b 



J b a (y) (3) 



where the current 



J b 1 (x) = i[A(x),d x A(x) 



q\{x)q a {x). 



(4) 



m is the 'current' quark mass and g is a coupling constant with the dimensions of mass. There 
is no term of the form q^Aq when non-dynamical fields are eliminated. This is essentially 
the dimensional reduction from 3 dimensions, except for the absence of terms of the form 
q^A-^-Aq and q^A-^-q. These theories have no ultraviolet divergences, though m undergoes a 

Ox Ox 

1 This dimensional reduction may well provide a first approximation to the gluon distribution of the proton 
measured in the Bjorken limit of deep inelastic scattering. 

2 The parallel transport operators between these points are trivial in a gauge where the null component of the 
the gauge field is set to zero. 
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finite renormalization (to 2 — > m 2 — g 2 N/ir). Heuristically, g 2 is related to the dimensionless 4d 
coupling constant g 2 via a factor of the transverse area of hadrons, g 2 ~ ^qcd5 , 4 • 

Before proceeding further, let us mention some literature on these theories. The adjoint scalar 
glueball spectrum was studied by numerical diagonalization of the hamiltonian by Klebanov and 
Dalley 4 . There is a large literature on discretized light-cone and transverse lattice approaches 
for which we refer to the review by Burkardt and Dalley Among the more analytical 
approaches, Rajeev, Turgut and Lee [H] studied the Poisson and Lie algebras of Wilson loops, 
open and closed string observables in the large- N limit of the hamiltonian framework. On the 
other hand, Rajeev formulated 2 dimensional QCD in the large- N limit (without the adjoint 
scalars) as a non-linear classical theory of gauge-invariant quark bilinears [7J. By developing 
approximation methods to solve this non-linear classical theory it was possible to analytically 
determine the quark structure of the baryon predicted by 2d QCD. This was shown to provide 
a first approximation to the non-perturbative quark structure functions of baryons measured in 
deep inelastic scattering [HI E] . 

Despite these developments, determining the expectation values of glue-ring and meson ob- 
servables analytically remains quite hard for two reasons. On the one hand, there are an infinite 
number of observables with ever increasing string length. Their dynamics is inter- linked via the 
factorized Schwinger-Dyson or loop equations. On the other hand, we are interested not just 
in the vacuum correlations but those in the ground state of the baryon, which has itself to be 
determined dynamically. 

We propose to use the approximate ground-state of the baryon determined in the large- N 
limit of 2d QCD as a starting point for the harder problem of determining the gluon correlations 
in the baryon. We have shown [HJ |S] that 

|*>= / dx 1 ---dx N e a '- aN ip( y x 1 ,---x N )q* ai ---q* aN (x N ) (5) 

for a factorized wave function (see Sec. |2|) provides a good approximation to the ground state 
of the baryon in the chiral and large- N limits. In a sense, the state |^) contains no gluons. 
Nevertheless, the presence of a baryon \^f) made of ./V quarks may be expected to deform the 
vacuum of the gluon field. The idea is to approximately determine the large- N expectation 
values of glue-ring observables Q in the background of this N quark baryon. Though \^f) is 
not the true baryon ground state of the action © , it is still an approximate ground state of the 
baryon number one sector of 2d QCD in the large- N limit. As such, it furnishes a variational 
approximation to the ground state of (j3J). In this manner, we hope to isolate a portion of this 
difficult problem that we have some chance of addressing analytically. 

In this paper we will take a very modest step towards the above-mentioned goal by studying 
a much simplified version via the path integral approach in the large- AT limit. To focus on 
the interesting baryon-part of the problem, we study a caricature of the action (j3J) without 
derivatives or non-local interactions and where fields are assumed time independent: 



S(A,q,q*) = tr J dx 



aA (x) — iq*(x)(m + aA (x))q(x) 



(6) 



Here a is the 't Hooft coupling with dimensions of mass, held fixed as N —* oo . Though an 
oversimplification, it allows us to isolate the dynamics associated to the degrees of freedom at 
the locations of the N quarks making up the baryon. After integrating out the quarks we get 
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an TV matrix model with an action involving not just traces of powers of the adjoint scalar, but 
also its inverse and determinant. To get a computable toy-model that preserves the essential 
picture of gluon correlations in a baryon background, we assume that the adjoint scalar field is 
equal at the positions of the different quarks. This leads us to a 1 -matrix model whose action 
is a quartic polynomial with a logarithmic modification. 



S(A) = tr 



A 4 - log [v + A 1 



(7) 



where the dimensionless parameter v is a ratio of quark mass to coupling constant. The loga- 
rithmic term encapsulates the effect of the A-quark baryon. Due to the significant truncations 
we make, this matrix model is not an approximation to the 1 + 1 dimensional field theory. It 
is the simplest toy-model where we may study the competing effects of gluon self-interactions 
and their interaction with quarks in a baryon. Finite-dimensional matrix models are often of 
independent interest. For instance, a slightly different logarithmic action log(l — A) + A appears 
in the Penner matrix model, which found applications to c = 1 string theory |lUj . 

Summary of Results: We solve our 1 -matrix model in the large- N limit by an extension 
of the methods used by Brezin et. al. jllj for polynomial potentials. This is possible because 
the derivative of the action is a rational function. We determine the free energy and glue-ring 
expectation values (moments of the adjoint scalar) in this caricature of a baryon. As v decreases 
below v c = 0.27, the 1-cut solution 3 of the matrix model makes a transition to a 2-cut solution. 
The reason is that the logarithmic term in the action (due to the presence of the baryon) makes 
small eigenvalues for the adjoint scalar energetically costly. The 1-cut solution could have been 
obtained by summing planar diagrams, but the 2-cut solution is analytically unrelated and could 
not be obtained that way. At the critical point v c , the leading discontinuity in free energy is in its 
3 rd derivative, while the 2 nd derivatives of the two- and four-point correlations are discontinuous. 
This large- A phase transition bears some resemblance to the Gross- Witten transition in the 
unitary one-plaquette model ° r that in the hermitian m 2 A 2 + gA A matrix model as m? / \J~g~ 
is made sufficiently negative All these phase transitions involve a jump discontinuity in 

the third derivative of free energy. One consequence is that analytic continuation from a weak- 
coupling treatment of gluons, starting from the vacuum {v = oo) would not be successful in 
predicting their behavior when strongly coupled to light quarks in a baryon (z/ — > + ). Moreover, 
moments of the adjoint scalar are enhanced in the baryon compared to the vacuum, since, in 
the two-cut phase, the support of the eigenvalue distribution excludes the origin. 

Organization of article: We begin with our ansatz for the baryon state in Sec. |2 and give 
a path integral formulation of the problem of finding glue-ring correlations in Sec. E3 Sec. E] 
describes the truncations and approximations that lead us to a matrix model for gluon correla- 
tions in a baryon. Sections El and El give the solution of the matrix model in the weak-coupling 
{y > v c ) and strong-coupling (y < v c ) phases respectively. The heavy quark limit, neighbor- 
hood of the critical point, and the chiral limit are treated in Seed Sec. |H1 presents a discussion 
of our results and open questions. 



3 One and two 'cuts' refer to the number of disjoint intervals of the real line on which the distribution of 
eigenvalues of A is supported. 
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2 Ansatz for baryon state 



A lesson from our study of 2d QCD in the large- N limit is that the ground state of the baryon 
is well approximated by a state containing N 'valence' quarks. Such a state may be built out of 

\a\;---;a%) = ^ (xi)gi 2 (x 2 ) • • • ql N (x N ) |0) (8) 

where q, satisfy canonical anti-commutation relations with respect to the Dirac vacuum |0). 
Ignoring flavor and spin quantum numbers, a general N -quark state is a linear combination 

\x£>) = £|o) = J dx 1 '''dx N 9(a\]---;aN)^ 1 {xi)^{x 2 )---^ N (x N ) |0>. (9) 

To represent a baryon, \I/ must be totally anti-symmetric in color 

*(a\;---;a%)=e a i- aN 4;(x 1 ,---,x N ), (10) 

and therefore represents a fermion if ip(xx, ■ ■ ■ , xn) is a symmetric function. For a path integral 
formulation, associate the grassmann-valued field q a (x) and its complex conjugate q* a {x) to the 
operators q a (x) and q\(x). Then the baryon state |^} becomes 

I*) -> jdx x ---dx N e ai - a " i>{x x ---x N ) q* ai {xi)q* a2 (x 2 )---q* aN {x N ). (11) 

By solving 2d QCD in the large- N limit we determined that ip{x\ ■ ■ ■ xn) is well-approximated 
by an iV-fold product of single particle wave-functions ip(x). A way to understand this is that 
once the antisymmetry in color is accounted for, quarks behave like N bosons which condense to 
the same one-particle ground-state. A good variational estimate turns out to be ip(x) = [fz^] 

which is the Fourier transform of the 'valence quark wave-function' ip{p) = 2y / 7re~ p 9(p > 0) . 
Though presented in the language of quantum many-body theory, the above ansatz may be 
obtained as the first term in a systematic field theoretic approximation method for the ground 
state of 2d QCD (see for details). One goal is to have a similar understanding of the gluon 
content of the baryon from the adjoint scalar field theory (J3J), which may then also be interpreted 
in terms of a relativistic many-body problem. 



3 Gluon correlations in the baryon state 

For any U(N) -invariant operator O , such as the glue-ring (^Q), the relation between expectation 
values of Heisenberg field operators and functional integrals in Euclidean space-time is 

(Z\0(A,q,q*)\® = (0\iom0) = ^JvAV q Vq*e- NS ^kO(A, q , q *)C 

where Z = J VAVqVq*e~ s{A > q ' q * ) tC (12) 

Here creates a baryon in the infinite past — T — ► — oo while £ annihilates it in the infinite 
future T . Using periodic boundary conditions in time we can assume that these operators are 
adjoints of each other and evaluated at a common time T. Specializing to our baryon state 
|f) = |^) given in (fTT|) . 

{%0{A,q,q*)\9) = I J ' VAVqVq*e- NS ^ A ^ J dx L ■ ■ ■ dx N e ai ... aN i>( Xl , ■ ■ ■ , x N ) 
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Xq a >(x u T)---q a »(x N ,T) 0(A,q,q*) J dyt • • • dy N e b ^ 

x i>* (yi , ■ ■ ■ , vn) <?& x (yi , T ) ■ ■ ■ 4 N (vn, t) ( 13) 

The functional integrals over Grassmann variables would vanish if xi ^ yi, so we take them to 
be equal and get 



(tf|0(A<Z, <?*)|*> = | j ' d Xl ---dx N \^{x lr --x N )\ 2 j 'vAVqVq*e- 

x q^( Xl ,T)---q aN (x N ,T) 0(A,q,q*) q* bl ( Xl ,T)---ql N (x N ,T) (14) 



y * p -NS(A,q,q*) ... e &1 "' 6j> 



where 

Z = J dx\ ■ ■ ■ dxN\ip(xx, ■ ■ ■ xn)\ 2 J T>AT>qT>q*e~ 

X ^ (x\,T) ■ ■ ■ q a « (x N , T) q* bl (x\,T) ■ ■ ■ q* bN (x N , T) 

dx\ ■ ■ ■ dxN \ip(x\, ■ ■ ■ xn)\ 2 Z(xi, ■ ■ ■ ,xn) (15) 



-NS(A,q,q*) b!-b N 
c ai---ajv c 



Now, we set aside the integral over x\, ■ ■ ■ ,xjy and focus on the functional integral 4 



Z{ Xl ,---x N ) = Jv[Aqq*]t 

q a ' (xi,T) ■ ■ ■ q a » (x N , T)q* bi (x\,T) ■ ■ ■ q* bN (x N ,T) (16) 



-NS(A,q,q*) M--b N 
c ai"-ajv c 



4 Matrix model caricature of baryon 

Physically, we are primarily interested in the shape and energy of the proton in its stationary 
ground state, which is probably best addressed in a hamiltonian framework. To mimic this in 
our path-integral approach, we will simply assume that fields are independent of time and ignore 
the dt terms in ©, which contain the symplectic structure in a hamiltonian approach. This 
can be regarded either as a crude passage to a kind of hamiltonian or as a further dimensional 
reduction. Moreover, it is clear that the locations of the quarks x±, ■ ■ ■ , xm play a special role in 
the path integral (fTT))) . We would like to isolate the contribution of those points. Of course, this 
is not possible with the Lagrangian (j3J), since it involves spatial derivatives as well as non-local 
interactions. However, in keeping with the modest goals of this paper, we work with an 'action' 
that is motivated by ©, but does not involve any derivatives or non-local interactions 

S(A,q,q*)=tr j ' dx\v{A) - iq* (m + aA 2 )q^ (17) 

In this caricature, V(A) = aA 4 stands in place of the current-current self-interaction of the ad- 
joint scalars in (jHJ). q*A 2 q models the quark-gluon cross term in the current-current interaction. 
We ignore the quartic term in quarks, which does not involve gluons. The primary effect of this 
term has already been taken into account in determining the baryon state 1^) (|11|) by solving 
2d QCD in the large- N limit 0. Thus (JTHJ) becomes 

Z(x u ---,x N ) = JvAVq*Vqe- N f dx[trV(A) - iq * {m+aA2)q] 

4 One complication that we do not address in this paper is that in calculating expectation values, the functional 
integral is not separately normalized, but only after integration over the points xi, ■ ■ ■ , xm ■ 
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xe^- aN e bl ... bN q bN (x N ) • • • q 1 ^ (x 1 )q* 1 (x 1 ) • • • q* aN (x N ) (18) 
The partition function can be written as: 

Z( Xl , ■ ■ ■ , x N ) = J VAe~ N tr / dxV (A) Zq (19) 
where Z q is the 'quark part' of the partition function 
Z q = J P^Pge^W^+^X^ 

= II / dq*(x)dq{x)e iNAq,{m+aA2)q e ai - aN e bl ... bN 

X 

q b »(x N )---q b i(x 1 )q* ai (x 1 )---q* aN (x N ). (20) 

In the second line we have discretized J dx to a sum J2 X A , where A ~ dx . Moreover, Z q may 
be factored into a 'vacuum part' and a 'baryon part' Z q = Z qv Z qb . 



Z qv = Y[ J dq*{x)dq{x)e iNAq ^ m+aA ^ q = [] det[NA{m + aA 2 (x))] 

X-f— X ^ X -^—X^ 

Z qb = e 6l ... bjv J] / ^*(^0d9(^)e iW(m+aA2), ? 6, (^)?: j (^)- (21) 



Now we use the result / dq*dqe' iq * Bq q b q* = [B det [B] to get 

{N f } %-N = z ai "' aN ^-b N d^B(x 1 )[B- 1 (x 1 )t 1 • • • detB(x N )[B-\x N )] b a \ (22) 

where B(xi) = m + aA 2 (xi) is an N x N matrix for each i = 1, ■ ■ ■ N. We can ignore the 
constant overall factor involving A . Thus Z qb depends on the gluon field at the location of the 
quarks. We factor the integral over A into a baryon contribution and vacuum contribution. 



f VAe -Ntr Jd X V(A) = -Q J dA{x)e -NAtrV(A) 

r N « 

x f[dA( Xi )e- NAt ^i=^ A ^\ (23) 
•* i=i 



Combining the quark and gluon vacuum contributions, we define the vacuum part of the partition 
function as 

Zvac(xi, ■ ■ ■ ,x N ) = J] J dA(x)e~ NA tr v ^ det [NA(m + aA\x i ))}. (24) 

X^-X-i 

Z vac is, for our purposes, the uninteresting part, so we consider the quotient 

N 



Zjxi ■ --xn) 

Z V ac{xi, • • • , Xn) 



f f[dA( Xi )e- Ntr ^ Av ^h a ^e bl ... bN 
J i=i 

xdetB(x 1 )[B- 1 (x 1 )] b a \---detB(x N )[B- 1 (x N )] b a 1 i (25) 

Here B(xi) = m + aA 2 (xi). This is an N -matrix model involving the N x N matrices A{x{), 
i = 1, ■ ■ ■ N . We do not yet have a way of approximately solving this multi-matrix model, so we 
will make some further assumptions that allow us to reduce it to a one-matrix model. 
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4.1 Reduction to a 1 -matrix model 



Suppose we assume that the adjoint scalar gluon field is equal at the positions of the N quarks, 
A(xi) = A for i = 1, ■ • • N . Then B(x\) = B{x2) = ■ ■ ■ = B(xn) = B and we can simplify Z q \,\ 



Then, 



Z qb = (detB) N e ai ... aN e b -- b -[B-X ■ • ■ [^C = (detB)^" 1 . (26) 
Z(xi ■ --xn) 



Z va c(xi, • • • , Xn) 



f dAe- NL ^ v ^ A \detB) N - 1 



— N 

dAe 



L tr V (A)- log det [m+ai 2 ] 



(27) 



where we have replaced N — 1 by N in anticipation of the large- N limit. iVA = L is a 
'length' of the baryon, assumed to have a finite limit as N — » oo and A — > 0. Recalling that 
= aA 4 , the action becomes 5(A) = tr [LaA 4 -log(m+aA 2 )} . Re-scaling A -» ^/(qL) 1 / 4 
we find there is only one dimensionless parameter v = myL/y/a, on which the observables of 
this matrix model depend non-trivially. The appearance of L is an artifact of our truncation 
and could probably be avoided. In the field theory Q, the relevant parameter would be m/g. 
Thus we arrive at the one-matrix model 

Z = J dAe~ N tr [^ 4 -i°s[^ 2 ]] . (28) 

v as a dimensionless ratio of quark mass to coupling constant. The absence of a quadratic term 
in A may be traced to the absence of a gluon mass term in (jSJ). Since A is hermitian, A 2 is 
positive. Thus, there is no difficulty in defining loghv + A 2 ] for v > 0. The expectation values 
of tr A n in this matrix model are what remain of the glue-ring expectation values. 



5 One-cut solution of quartic + log matrix model 

The rest of this paper will be devoted to a study of the large- iV limit of the one- matrix model 
with action and free energy 



S(A) = tr 



A 4 - log(i/ + A 2 



, E(v) = - lim -^-log [dAe- NS ( A \ (29) 

JV^oo N 2 J 

Glue-ring expectation values 

G n (u) = lim ( [ dAe- NS W±-tvA n )/ [ dAe~ NS ^ (30) 

N^co \J N J J 

are given by moments G n = J p(x)x n dx of the eigenvalue density. p(x) must minimize the free 
energy 

E[p] = J S{x)p{x)dx — J dxdyp(x)p(y) log \x — y\; E(u)=mxn E[p] (31) 

where S(x) = x 4 — log (u + x 2 ) . (a) For v sufficiently large, Six) is convex from below and we 
expect p to be supported on a single interval, (b) For v sufficiently small (but positive), S(x) 
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is shaped like a Mexican hat (see Fig. H3) and we expect p to be supported on a pair of intervals 
located near the minima of S(x). These two cases are treated in Sec. Eland 13 For p supported 
on a single interval, the Mehta-Dyson linear integral equation for an extremum of E\p\ is 



S'(x) = Ax 3 



2x 



V + X 1 



2V 



2a 



p(y) 



2a X - y 



dy 



2a < x < 2a, a > 0. 



(32) 



p is subject to positivity p(x) > and normalization J^ a p(y)dy = 1 conditions. Since S(x 
is even, p must be even. It is convenient to introduce the generating function of moments 



F(z) 



(33) 



n=0 



I -2a Z-y 

(a) F{z) is analytic on C\[— 2a, 2a] . (b) F{z) ~ ^ as \z\ — > oo which follows from normalization, 
(c) F(z) is real for real z outside [—2a, 2a] . (d) When z approaches [—2a, 2a], using (|32]> . 



1 2x 

Fix ± ie) = -(4x 3 =)^mp(x) =^ p(x) 

2 v + x z 



1 

2vrf 



(F(x - ie) - F{x + ie)) (34) 



(e) F{z) is an odd function of z since p is even. By analogy with the case of polynomial S(x) 
we expect that if there is an F{z) satisfying these conditions, it is unique. Existence, however, 
is not guaranteed. Indeed, we do not expect a 1-cut solution for sufficiently small v > 0. In a 
region of validity which is to be determined, the following ansatz for F(z) is consistent with the 
above requirements 



F(z) = 2z 3 



v + z- 



+ R(z)Vz 2 -4a 2 . 



(35) 



R(z) is a rational function to be chosen so that these conditions are satisfied. In particular, we 
must pick R{z) to cancel the poles at z = coming from . So R(z) = > where 

the polynomial P{z) does not have zeros at z = ±i\fv . Moreover, to cancel the linear and cubic 
terms in S' (z) as \z\ — » oo, we need to pick P(z) to be a quartic even polynomial. Thus 



F(z) = 2z 3 



V + z 2 



+ 



2 i 4 

a + "fz + ez 
(u + z 2 ) 



V z 2 - 4a 2 . 



(36) 



We need to determine the parameters a, a, 7, e. Though it may appear that the third term is 
even in z, it is actually odd due to the square root. As z — > 00, 



F(z) 



(2 + e)z 3 + 
+ 



z z* 



7 - e(2a 2 + v 
1 + a - 2a 4 e - + eis 2 - 2a 2 (j - ev) 

Requiring F(z) ~ i fixes 0,7 and e 

e = -2; 7 = -(2z^ + 4a 2 ); a = 2 - 12a 4 - 4^a 2 . 
To fix a we must require analyticity of F{z) as z 2 — > — ^. This is ensured if 

- 2 + (a + 7Z 2 + ez 4 )^^ 2 - 4a 2 



(37) 



(38) 



(39) 
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vanishes as z — > ±«i/z/. In other words, 



± i\J~u + (a — 71/ + ev )^J-v-\a 2 = (40) 

or —1/ + (—2 + 12a 4 ) 2 (4a 2 + ^) = 0, which is a quintic equation for s = a 2 , 

576s 5 + 144i/s 4 - 192s 3 - 48z^s 2 + 16s + 3u = 0. (41) 

We are guaranteed at least one real solution for a 2 . The physically allowed solutions are those 
with a 2 > 0. For v > 0, we find that there are two solutions a 2 > only the smaller of which 
leads to p(x) normalized to 1 (the other has p normalized to 3). For this value of a, the density 
of eigenvalues is 

p( x ) = -(-) R(x) V±a 2 - x 2 = - ( - V Q ~, + J X ) y/4a 2 - x 2 (42) 
\ttJ \ttJ (v + x z ) 

with a, 7,6 given in (|38[). 



5.1 Moments 

The analogue of glue-ring expectation values (|30|) can be obtained from the Laurent series for 
the moment generating function F(z) 

F(z) = 2z 3 + — jl — Vz 2 - 4a 2 = V — ^ (43) 

w i/ + z 2 (z^ + z 2 ) v ; 

v ' 71=0 

The odd moments vanish G2 n +i = and the even ones are 
G 2 = 40a 6 + I2va 4 - 4a 2 - u, 

G 4 = 60a 8 - 24ua 6 - 4(1 + 3z^ 2 )a 4 + 4i^a 2 + u 2 , etc. (44) 
where a{y) = yfs is the physical solution of 1)41(1 . G^v) an d ^4(1/) are plotted in Fig. ^andEl 



5.2 Free energy 

Let us define the entropy ^1] as 

x = v 



p(x)p(y) log \ x — y\dxdy. 



(45) 



The large- N free energy is the Legendre transform of entropy, E = J dxp(x)S(x) — x- An 
expression for E involving only single integrals can be obtained using the Mehta-Dyson equation 



\s'{z) = V [ -^-dy; z £ supp(p). 
2 J z-y 



p{y) 



y 

Integrating with respect to z from xq to x, both of which lie in the support of p, 



S(x) - S(x ) 



V J dy p(y) 



log \x-y\- log \x - y\ 



(46) 



(47) 
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Two-Point Function G2 (v) 




Figure 1: Two-point glue-ring expectation value G%{y) plotted against logzA The asymptotic 
behavior for heavy quarks or weak-coupling {v — > 00, 1-cut phase) is separated from the behavior 
for gluons strongly coupled to light quarks ( v — ► , 2-cut phase) by a third-order phase transition 
at v c . G2 grows as we go from gluons in the vacuum (u = 00) to those in a baryon strongly 
coupled to light quarks (small v) since p{x) transforms from uni-modal to bi-modal. 



Now multiply by p{x) and integrate with respect to x, 
1 



- J dx p(x) 



S(x) - S(x Q ) 



V 



dxdyp(x)p{y) 



log \x-y\- log \x - y\ 



(48) 



Thus 



S(x) — log \x — xq\ 



E = -S(x Q ) + V J dx p{x) 

1 1 f 2a 

E{u) = S(x Q ) + -G A -V dxp{x)\og\{v + x 2 ){x 2 -xl)\. 

2 2 Jo 



(49) 



xq is arbitrary provided it lies in the support of p. The xq independence of E follows from the 
Mehta-Dyson equation and was also verified numerically, xq = is most convenient for us, so 



1 1 /" 2a 

E(y) = — log v + -G A (v) -VI p(x) log {vx 2 + x 4 )dx. 

2 2 Jo 



(50) 



The new ingredient in free energy not determined by polynomial moments G n is a sort of 
logarithmic moment. We could not evaluate it in terms of known functions for the 1-cut p{x) 
but the integration is easily performed numerically and plotted in Fig. [21 



5.3 Domain of validity of one-cut solution 

For which v > is the 1-cut solution valid? So far we have not imposed the p > condition. 
The 1-cut solution 

p(x) = -(1/tt) , P ^ X \, Via 2 - x 2 ; P(x) = a + 7 x 2 + ex 4 (51) 
(y + x z ) 

will break down if p < 0. The boundary of the region of validity is given by v for which 
p(x) = for some |x| < 2a. The most obvious way in which the 1-cut solution breaks down is 
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Free Energy E(v) 
1.75(- 




Figure 2: Free energy E{y) versus logz/. The 2-cut phase lies to the left while the 1-cut phase 
lies to the right of the critical point at v c . For heavy quarks or weak-coupling (y — ► oo), E{y) — ► 
— log v — .996 as shown in Sec. l7.1l and illustrated by the linear asymptotic behavior to the right 
in the plot. In the chiral limit (Sec. 17.3]) . free energy approaches a constant E(0) ~ 1.751. 

when p(0) = and we transit to a 2-cut solution without any support in an interval containing 
x = 0. p(0) = = implies a = 0, since a > 0. a = 6a 4 + 2va 2 -1 = 0. Regarded 

as a quadratic in a 2 , the unique positive solution is 



a 2 = -^ + ^±^>0 for u>0. (52) 
6 6 

This is the critical value of a(v) at which the 1-cut solution vanishes at the origin. Substituting 
a 2 into the quintic (14111 we find a condition on the transition point v = v c 



27 - 8z/ + 48i/ J(6 + v 2 ) + 8z/ 2 (-9 + uJ{6 + z^ 2 )) 



0. (53) 



Assuming v > 0, the only possibility is for the factor in parentheses to vanish. Upon simplifi- 
cation it becomes a quadratic equation 48z^ 4 + 368i^ 2 — 27 = 0. The positive solution is 



/13^/l3-46 . , 

Vc = i^~Y2 PS0 ' 27 - (54) 

Are there any other transition points, i.e. does p(x) become negative for any < |x| < 2a? 

For example we find that there is no value of v > for which P(2a) = 0. Indeed, the condition 
P(2a) = implies that 30a 4 + 6isa 2 - 1 = or 

2 v V9v 2 + 30 , . 

a= . (55) 

10 30 y ' 

However, when this is substituted into the quintic (|41l) . the condition 

3200 



4875z^ - 600^ 3 + 72^ 5 = 24^ 4 240v 2 V30 + 9v 2 (56) 



has no solution for v > 0. Based on the shape of S(x) we expect the 1-cut solution to be valid 
for all v > u c , and we have checked that this is indeed the case. 
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6 Two-cut solution of quartic + log matrix model 



For small is, S(x) = x 4 — log[f + x 2 ] develops a repulsive core near x = 0. For v < u c , we 
expect the 1-cut solution of the Mehta-Dyson equation 



S'(x) = 4x 3 ^ =2T ( -?QLdy, x e supp(p) (57) 

v + x z J x — y 

to make a transition to a 2-cut solution supported on [—2a,— 26] U [26,2a]. The generating 
function of moments F(z) = J jr^jdy enjoys the same properties (a) -(e) as before with [—2a, 2a] 
replaced with [—2a, —26] U [26,2a] where a > 6 > 0. An appropriate ansatz for F(z) is 

F(z) = ^S'(z) + R(z)^/(z 2 - 4a 2 )y/(z 2 - 46 2 ) 

= ^^^^^ <58> 

(3,5,a,b have to be fixed using analyticity and asymptotic behavior of F(z). As \z\ — » oo, 

-► ^2 + ^ z 3 + (j3 - 5{2a 2 + 26 2 + v)^ z 

+ ( - 1 - 2/3(a 2 + 6 2 ) - 25(a 2 - 6 2 ) 2 - v{$ - 25{a 2 + 6 2 )) + Su 2 ^j 1 + O(-I) (59) 

The requirement F(z) ~ j + 0(l/.z 3 ) implies 

6 = -2; /? = -2(2a 2 + 26 2 + i/) and 6(a 4 + 6 4 ) + 4a 2 6 2 + 2v(a 2 + b 2 ) - 1 = 0. (60) 

The condition that F(z) be analytic at z = ±iyfv implies {(3 — v5)^J{v + 4a 2 ){u + 46 2 ) + 1 = 0. 
Substituting for (3 and (5, we are left with a pair of algebraic equations for a and 6. 

6(a 2 + 6 2 ) 2 - 8a 2 6 2 + 2v(a 2 + 6 2 ) - 1 = 
4(a 2 + 6 2 )y / 16a 2 6 2 + 4^(a 2 + 6 2 ) + v 2 - 1 = 0. (61) 

Let s = a 2 + 6 2 and p = a 2 6 2 be the sum and product, then 



6s 2 - 8p + 2vs - 1 = and 4sy / 16p + 4vs + v 2 -1 = 0. (62) 
We can eliminate 

1 



P= 4 



2 1 

3s + vs 



(63) 



and get an algebraic equation for s, 4s^l2s 2 + 8vs + (y 2 — 2) = 1, which has at most one 
positive solution s for v > . Squaring it we get a quartic equation 

12s 4 + 8z/s 3 + {u 2 - 2)s 2 - — = 0. (64) 

This can be solved in lengthy but closed form and the unique positive s selected. From s, we 
get p as well and 



2a = ^2s + 2y^ 2 -4p, 26 = \j 2s - 2<J~^ - 4p. 



(65) 
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Four-Point Function G4 (v) 




Figure 3: Four-point glue-ring expectation value G±{v) versus log^. The 1-cut solution for 
v > v c and 2-cut solution for v < v c are not analytic continuations of each other despite 
appearances. The asymptotic values at v = 0, 00 and at v c shown on the vertical axis are 
obtained analytically in Sec. [Tj 

The eigenvalue density, supported on 26 < \x\ <2a is 
p[x) = -^M,/(4a 2 -x 2 )(x 2 -462) = + ( 4 « 2 + ^ + 2 ") x ) / (4fl2 _ x 2 )(x 2 _ 4b 2), (66) 

7T v 7r(z/ + X 1 ) v 

Using the Laurent expansion of F(z) , we get the moments G2 n +i = 0, 
G 2 = 16s 3 - 64ps - 8vs 2 - Av 2 s + u, 

G A = 64p 2 + 36s 4 - 32ps(5s - z;) - 8s 3 i/ - z; 2 + 8s V + 4sz/\ ••• (67) 

s and p can be eliminated using the solution of the quartic equation (|64j) . The results are 
plotted in Fig. [I] and [U 

6.1 Free Energy 

The free energy may be expressed in terms of single integrals 

E{v) = ^S(x ) + \g± - l -V j dxp(x)\og\{u + x 2 ){x 2 -x 2 )\ (68) 

where xq is any point in the support of p. To obtain (|68[) we need to pay attention to the fact 
that the Mehta-Dyson equation (|57|) is valid in two disjoint intervals. Begin by integrating the 
Mehta-Dyson equation with respect to z from x + to x where x + ,x S [26,2a] to get 5 

~[S(x) - S(x + )] = J dyp(y)[log\x - y\ - log\x + - y\] for x,x+e[26,2a] (69) 

Multiplying by p(x) and integrating with respect to x from 26 to 2a and simplifying gives 
f 2a f 1 f 2a <S(x+) If 

j2b dx J dyp ^ p ^ log \ x ~y\ = 2 J 2b dx p( x ") s ( x ) — 4 — h 2 J dyp ^ log \ x+ ~ y ' ^ 70 ^ 

5 The principle value prescription is implied where necessary and for brevity will not be explicitly indicated. 
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X 



for x+ £ [26, 2a] . When the limits of integration are not specified, the integral is over supp(p) . 
Similarly, for x_ £ [—2a, —26] we get 

/■-2b r l r-2b S(X-) 1 f 

J 2 dx J dyp ^ p ^ log ' x ~ y ' = 2 7 2 cfa/9 ( x ) 5 '( 2; ) — 4 — h 2 7 d yp(y) lo % \ x - ~ y\- ( 71 ) 

Adding these two, we get for the entropy (d5)> , 

X = i / dxp(z)S(x) - - + i y dypfo) log |(x + - y)(x_ - y)\ (72) 

where x-t are in the positive and negative part of supp (p) . By choosing x_ = — x + we simplify 
matters using the fact that S(x) is an even function 

If 1 1 f 

- I p(x)S(x)dx - -S(x + ) + - dx p(x) log |x 2 - x\\. (73) 

Now observe that in this expression, the sign of x+ does not matter, so we can call xq = x+ and 
pick it anywhere in supp(p) . Recalling that E = J p{x)S{x)dx — x we get the advertised expres- 
sion 1)68(1 for the 2-cut free energy. Though the complexity in evaluating E(v) has been reduced, 
we have not been able to find the above logarithmic moment in terms of known functions. The 
numerically evaluated free energy E(v) is plotted in Fig. [21 

6.2 Phase transition to one-cut solution 

We expect the 2-cut solution to make a transition to the 1-cut solution when the intervals 
[—2a, —26], [26, 2a] merge, i.e. 6 = 0, which implies p = and s = a 2 . Inserting in (|63|) gives a 
quadratic equation 3s 2 + vs c — \ = whose positive solution s c = a 2 . = — | + ^ U ^ +G is the same 
as the value of a 2 at which the 1-cut solution breaks down. When s c is substituted in (|64l) . we 
get the same phase transition point v = v c as before (|54[l. The 2-cut solution for v < v c takes 
over when the 1-cut solution breaks down. 



7 Special cases: weak-coupling, critical point and chiral limit 

7.1 Heavy quark or weak-coupling limit v — * oo 

In our toy-model, when the quarks are very heavy or the coupling constant is small, the self- 
interactions of the gluons dominates their interactions with the quarks in the baryon. This is 
the limit v — > oo where the action S(A) — > tr A 4 up to an additive constant. It is as if the 
gluons do not feel the presence of the baryon and we return to calculating vacuum correlations 
of glue-ring observables. This limit lies in the deep end of the 1-cut phase, where calculations 
simplify. The quintic equation (|41|) for the limits of supp(p) reduces to a quadratic equation 
48s 4 — 16s 2 + 1 = whose physical solution is s = \/3/6. The limiting eigenvalue density is 

1 



p(x, v ->■ oo) = -(4s + 2x 2 ) V4s - x 2 , \x\<2y/s, (74) 



7T 



The odd moments vanish while the even moments and free energy (|49|) are 

2 n+1 (n + l)r(n + i) 2V3 1 

3 n /Virr(n + 3) 9 4' 



14 




.1 .2 v c .3 .4 .5 

Figure 4: Derivative of the two-point correlation G<i(v) has a kink at v c = 0.27 indicating that 
its second derivative is discontinuous. 



Second Derivative of G4 (v) 



1.8 
1 . 6 
1 . 4 
1 .2 

1 

. 8 




.1 .2 v c .3 .4 .5 
Figure 5: Jump discontinuity in the second derivative of G±{y) at the critical point v c . 

11/" 1 

£ = -log^+-G 4 / dx p{x) log x 2 = - log v+ -(3 + log 144) k, -log v + .996.(75) 

11) 8 

These limiting values are seen to agree with the numerically obtained behavior of E{y) , G^iy) 
and G&{y) plotted in the Fig. E an( i El f° r a wide range of values of v . 



7.2 Neighborhood of the phase transition 

At u c , the 1-cut solution, valid for high quark mass m or weak-coupling a makes a phase 
transition to a 2-cut solution, which is valid for low quark mass or strong-coupling {y < u c ). 
In the immediate vicinity of the phase transition, observables are more easily evaluated than 
generically. We find that the eigenvalue density is continuous across the transition, the critical 
eigenvalue density is (see Fig. |S} 



Pc{x) 



{(7V13 - 11)v c x l + (13V13 - 46)x 4 } 

187Tl/2(zy c + x 2 ) 



13 + 20> c - 9x 2 



(76) 



and is supported on [— 1a c ,1a c ] where a 2 c = y ( ^ 2 ) , a c ~ .605. Thus, all moments and the 
free energy are also continuous. The next question is whether their derivatives are discontinuous 
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p c (x) and S c (x) 



1 

\ 


1 




1 


\ 


. 8- 




1 




\ </■ 6 


- \ / 


A 


\ 


\/ . 4 




\ 


/\.2 


- / \ 


/ 



-1 -0.5 0.5 1 



Figure 6: Eigenvalue density at the critical point p c (x) , the continuous curve. The dashed curve 
is the critical action S(x,v c ) plus a constant chosen to make it vanish at its minima. While p 
is peaked around the minima of S , the spread allows it to maximize entropy % (E5J)- 



across the transition. We find that the second derivatives of G2 and G4 have jump discontinuities 
across the phase transition. We have calculated this behavior analytically. To find the derivatives 
of the correlations (|44j l .(|67 |) at v c we need the values and derivatives of a,s and p at v c . The 

critical values v c = V^v^s- Jts ,s c = a 2 ,p c = have already been determined without much 
trouble. For their derivatives, it is again not necessary to solve the quartic Q64JI and quintic 
(j41|) equations, but suffices to solve linear equations obtained by differentiating these at v c . For 
example, suppose we want G' 2 (z/r), which corresponds to the approach from the 2-cut phase. 
From (fo7|) 

G' 2 (v) = 48s 2 s' -64ps' - 64sp -4u 2 s - 8su - 8s 2 - 16ssW 1 (77) 
Differentiating (|64|) and solving the linear equation for s'(v) gives 

1 , s —vs — 4s 2 . 

S'M = n n- (78) 

v 7 -2 + v 2 + 12z^s + 24s 2 v ; 

p'{v) is determined similarly using ()63|) and evaluated at In this manner we get 

± (WTj + 19) _ , ± vg-5 ^ 
G 2(^ c ) = 27 ^ c - 66; G2 ^ c ) = 3 ~- 46 



_ 4 (205^13-122) _ 8 (1669VB-5858) 

2V c 7 9 V 39 2V c 7 51 V 39 v ; 



and 



(4706_U50^^ 844; (2522^14^ , 1J0B (m) 

bod 351 

These agree with the numerically determined correlations plotted in Fig. |I] and Fig. [5J 
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Second Derivative of Free Energy 




The free energy E(y) and its first two derivatives are continuous at v c . E" 1 (v) has a jump 
discontinuity at v c . For example, to calculate E'(v^) , we differentiate the integral representation 
for 1-cut free energy (jl9|) 

1 1 [ 2a 

E(y) = — log v H — G4 — / dx p v {x) log (yx 2 + x 4 ) 

2 2 Jo 

=* = -^ + ^)- rdx{^ + ^log(^x 2 + x 4 )j. (81) 

can be got from the explicit formula l|42|). We omitted the term involving the derivative of 
the upper limit of integration because p(2a) = 0. Evaluating at u c , using the above result for 
G'±(v) and doing the integral numerically gives us E'(i/^). Proceeding along these lines we get 
E(y c ) « 1.36, E'{y c ) w -1.32,E"(v c ) » .93,E"'(v-) » -.7, E"'{v+) « 7.05 This is illustrated in 
Fig. d We conclude that the phase transition is of third-order. 



7.3 Chiral limit or strong-coupling limit v — ► 

In the limit of massless quarks or strong-coupling, v — > + , which lies in the 2-cut phase. The 
action becomes S(A) = tr [A 4 — log(^4 2 )] and the gluons maximally feel the presence of quarks in 
the baryon. The solution again simplifies. The quartic equation for s (|64|) becomes a quadratic 

with solution s = \ \J 2+ ^ and p = "^" 2 . From this we get 

= (2x 2 + 4a) / 2 _ ^ 4 _ — 

7TT » 



and moments G2 = - — g — -y 4+ 3 ~ -794 and G4 = | which are shown in Fig. ^and|21 The 
free energy (jfi8j) has a finite limit E{y = 0) ~ 1.751 as illustrated in Fig. El 



8 Discussion 



We have studied a very simple matrix model for glue-ring correlations in a baryon in the limit 
of many colors. It was obtained as a caricature of the dimensional reduction of QCD to 1 + 1 
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dimensions. Our main finding is that there is a third-order phase transition that separates a 
phase where gluons are weakly coupled to heavy quarks ( v > u c ) from one where the quarks are 
light and strongly coupled to gluons (0 < v < v c ). u, a dimensionless ratio of quark mass to 
coupling constant is the only parameter of the model. While for v > u c we have a 1-cut solution 
of the matrix model, for v < u c we have a 2-cut solution. The gauge-invariant observables (glue- 
ring expectation values) are described in these two phases by two different analytic functions of 
v that disagree beyond their first derivatives at v c (See Fig. ^ El El EJ) ■ Moreover, the case of 
gluons in the vacuum i.e. where the baryon is absent, corresponds to v — > oo, which is deep 
inside the 1-cut phase. The physically interesting value of v most likely is small and lies in 
the 2-cut phase, since current quarks in the proton are very light compared to Aqcd • Thus, 
the vacuum correlations of gluons are likely to be separated by a phase transition from those 
in a baryon state. Moreover, from Fig. ^ and 01 we see that gluon correlations are enhanced 
inside the baryon compared to their values in the vacuum (y = oo) . This reflects the growth of 
moments of the distribution of eigenvalues, as they make a transition from being clustered about 
the origin to being supported on a pair of intervals excluding the origin. Though this is very 
far from explaining why about half the proton's momentum is contributed by gluons, it does 
indicate that the qualitative features of gluon correlations in the vacuum can be quite different 
from those in a baryon state. Even if we could use a weak-coupling expansion to describe a 
bound-state like a baryon, the phase transition would invalidate analytic continuation to the 
physically relevant baryon containing gluons strongly coupled to light quarks. 

The wider applicability of our results is called into question by our approximations and 
truncations. The sharp phase transition is an artifact of N = oo . For finite N , the matrix model 
has finitely many degrees of freedom and cannot display non-analytic behavior. Nevertheless, 
this is probably the most benign of our approximations. The finite N theory should display 
qualitative differences between the two regimes. Absence of space-time derivatives and non-local 
interactions due to longitudinal gluons are the more significant shortcomings of our model. This 
toy-model has given us a cartoon of how the theory may behave as v is varied. A hamiltonian 
approach may have a better chance at shedding light on the matrix field theory or matrix 
quantum mechanics version of this problem. We hope the proposal of treating gluons in a 'fixed 
baryon background' |\&) (|11|). along with other new ideas will help simplify the matrix field 
theory in order to better understand the emergent bound-state structure of gluons in a nucleon. 

Some questions for future work are collected here, (i) Can we use a variational or other 
approximation method to understand this phase transition? Such an approach has a better 
chance of generalizing to multi-matrix models, (ii) Besides glue-ring correlations, we are also 
interested in open string correlations in a baryon state. Can we get a zero-dimensional toy- 
model for these as well? (iii) Can we shed any light on the multi-matrix model that arises 
when we do not assume the transverse gluon field to be equal at the positions of the quarks? 
(iv) What is the relation between the gluon distribution function extracted from experimental 
data and the gauge-invariant glue-ring variables that would come from solving the matrix field 
theory? One suspects that the gluon distribution is essentially determined by a two-point glue- 
ring expectation value. What is the simplest truncated form of the matrix field theory where 
such a gluon distribution function can be estimated? (v) In [Tl], we obtained algebraic and 
probabilistic characterizations of the entropy x whose Legendre transform is the free energy of 
matrix models. However, the formula for x i n generic multi-matrix models is quite complicated. 
Is there any multi-matrix generalization of the trick (|49|) of using the Mehta-Dyson 'equation of 
motion' to reduce double integrals to single integrals? (vi) Can we find a variational principle 
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that determines the closed and open string observables of our matrix field theory or some finite 
dimensional truncation thereof? Such a variational principle for closed string observables in 
multi- matrix models was found in |14j . 
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